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Abstract 

Relaxor ferroelectrics, a typical example of which is PbMg 1 / 3 Nb 2 /303 (PMN), a per- 
ovskite in which the Mg 2+ and Nb 5+ randomly occupy the octahedrally coordinated 
site, have a very high dielectric constant for a very wide region in temperature and 
have no ferroelectric transition except in applied electric fields. They are technolog- 
ically important materials for transducers and ceramic capacitors in the electronic 
industries and also present an interesting physics problem. Cowley et al. in a 
detailed review of relaxor ferroelectrics have commented that although they were 
first synthesized over fifty years ago and their properties have been well explored, a 
satisfactory theory explaining their properties has not yet been formulated. We con- 
struct a minimal model for relaxor ferroelectrics based on Onsager's results on models 
with dipolar interactions, and the short-range displacive transition model with which 
pure ferroelectrics were understood, together with quenched disorder. We formu- 
late the simplest necessary approximate solution of the model with which each the 
macroscopic properties and recently measured correlation functions and their inter- 
relationship are understood in terms of a small number of parameters in the model. 
Specifically, we calculate and compare with experiments the structure factor as a 
function of temperature and quantify the various regimes in the dielectric constant 
of relaxors, the Curie- Weiss regime above the Burns temperature Tb, and the growth 
of unusual correlations between Tb and the temperature T max where the dielectric 
constant is maximum, recently discovered by neutron diffuse scattering. 
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The properties of the relaxors may be summarized as follows: in zero external electric 
field they have a region of collective dielectric fluctuations which is bounded at the upper 



end by a temperature Tg, the Burns temperature |2j, above which the susceptibility has 

n 

a Curie- Weiss form [3|. A broad region extending down to T = below Tg is marked by 



a temperature T max where the susceptibility is maximum and has a width in temperature 
AT max . Both Tmax and AT max depend on the frequency u at which the dielectric susceptibil- 
ity is measured 



4j. Very significantly, neutron scattering experiments js| have revealed that 
the structure factor S q has unusual power-law form in each regime of temperature, which 
is unlike that in the random field or the random bond disorder models with short-range 
interactions. Typically only such interaction models (which are sufficient to describe pure 
ferroelectrics) with quenched disorder [8| have been used to describe relaxor ferroelectrics. 
These models in three dimensions for weak disorder also have long-range order at a finite 
te,„pe r at ure fl, which ,s not observed in the c_l r e.axor PMN Q. 

In formulating the minimum necessary model for relaxor ferroelectrics, it is important to 



recall Onsager's result js] that unlike in the Clausius-Mossotti or Lorentz approximation, 
dipole interactions alone do not lead to ferroelectic order except at T = 0. Also that 
pure ferroelectric transitions were understood with the realization |7[ that they are soft 
transverse optic mode transitions due to short-range interactions that lead to structures 
without a center of symmetry so that ordered dipoles are induced. The first point is not in 
practice important for pure ferroelectrics which are well described by a mean-field theory 
for the dielectric constant, and is often a first order transition only below which the dipoles 
are produced. But as we show below, in relaxors the random location of defects acts in 
concert with the dipole interactions to extend the region of fluctuations to zero temperatures. 
Therefore dipole interactions must be a necessary part of the model. The essential physical 
points in the simplest necessary solution is to formulate a theory which considers thermal 
and quantum fluctuations at least at the level of the Onsager approximation and random 
field fluctuations at least at the level of a replica theory. 

We focus on the relevant transverse optic mode configuration coordinate Ui of the ions 
in the unit cell % along the polar axis (chosen to be the z-axis). itj experiences a local 
random field hi with probability P({hi}) due to the compositional disorder introduced by the 
different ionic radii and different valencies of Mg +2 , Nb +5 in PMN. The model Hamiltonian 
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IS 

H = E \wi + v{Ui) ] VijUiU > ~ E E ^ Ui - E ^ 

where is the momentum conjugate to ttj, M is an effective mass and we have included 
an infinitesimal external field Ef xt . We assume the h^s are independent random variables 
with zero mean and variance A 2 . V{uj) is an anharmonic potential, 

Vim) = \u\ + luj, (2) 

where k, 7 are positive constants. Vij is the dipole interaction, 

(3) 

0, Ri = Rj, 

where a is the lattice constant, z* is the effective charge and is the z-component of -Rj. 
The Fourier transform of the dipole interaction v q = ■Vije iq '( Ri ~ Rj > , is non-analytic for 
q — > 0. For future use, we denote Vq (= 47inz* 2 /3) for the q — > component of v q in the 
direction transverse to the polar axis (q _L z), n the number of unit cells per unit volume, 
and a the lattice constant. 

The Hamiltonian (TjQ) presents (i) long ranged (anisotropic) dipolar interactions; (ii) com- 
positional disorder; and (iii) anharmonicity. We now present the procedure used to find the 
correlation functions for (CQ). We use the quasi-harmonic approximation, similar to that used 
for pure displacive ferroelectrics to treat anharmonicity [9|. This reduces the problem to that 
of harmonic oscillators with a renormalized stiffness which is determined self-consistently. 
Disorder is treated by an approach which may be termed the "poor man's replica method" , 
and has been used to study two and three dimensional magnets with quenched random 



fields 10[. The Onsager approximation, which may be considered the leading fluctuation 



correction over the mean-field approximation 11] is used for the dipolar interactions as well 
as for the general statistical-mechanical treatment of the problem. 
In this approximation, 

H « E H * = E ( H " - h ^ - Er ity ui) , (4) 

i i 

with H9 = p 2 i/{2M) + V(ui) - (X/2)ul E- avity is the Onsager cavity field, 

ET mtV = E % («i>T " A <«*>t + E t Xt - (5) 



{Xi) T denotes the thermal average of any quantity Xf, similarly we will denote by ((Xj) T ) , 
the configurational average over the quenched random fields hi taken after the thermal 
average is taken. 

Onsager calculated the parameter A for the simpler problem of dipolar forces alone using 
continuum electrostatics. More generally, one may determine A by a self-consistent procedure 



such that the response functions obey the fluctuation dissipation theorem 12[. To find the 
self-consistent response functions, consider the linear response to the cavity field E^ amty 9]. 
Then 

( Mj > T = 0^)£r % , (6) 

where (phX 00 ) * s the dielectric susceptibility for the problem with Hamiltonian, Hf — U{h{. 
To take the average over configurations, we assume that the effects of the random fields hi 
decouple in <fihi(u) and (ui) T . For our model, as will be shown below, (</>hi{v)) c is independent 
of i and therefore has only a q = component denoted simply by (^(w)}. Then taking a 
Fourier transform, we get the formal expression for the susceptibility function, 

_ « Ug >T> C _ 1 

Next we determine <fih(u). The equation of motion for it, due to H — h(Ui and including a 
damping force characterized by T are 

Mil, = + + hi + ET\t) - Tin (8) 

We now linearize OH]) by considering fluctuations 5ui around the (static) expectation value 
of Ui, Ui = ) + Sui, to get that 

Mil' — uj z + tui 



where Mtt' 2 = k - A + 3 7 (((K)t) + 



The fluctuations ( (5uf)^J of the H — hiUi problem are easily shown to be, 

h ( (3hTl'\ I 2\ // \0\ 2 



^t)c = jm- coth {^r) + (<«<>*). - «"'>^ (10) 

We finally find the susceptibility (j7J), 



where Q q is the vibration frequency of Ui for the full problem, 

Mnl = M(pk) 2 +(y£-v q ) (12) 

and Qq is the q = component of the phonon frequency in the direction perpendicular to 
the polar axis (q J_ z), 

M (Q^) 2 = MQ' 2 - «-A). (13) 
We now determine the parameter A by enforcing the fluctuation dissipation theorem, 

(»;) c -(W: 2 ) c = £/-oth Imix^)} . (14) 

The summation over q extends over the first Brillouin zone. To close the system of equations 
the polarization (uj)^ must itself be related to the random fields through the susceptibil- 
ity ffTTT) . For a fixed realization of disorder 10|, 



(u i )° T = J2{xi j ) c h j , (15) 
j 

where {Xij{fy) c is the zero-frequency susceptibility averaged over compositional disorder with 
Fourier transform = Ylij (Xij(0)) c e iq ' ( - Ri ~ R ^ . Using (hj) c = and (hihj) c = A 2 ^-, 

««*>?■>„ = 0, ((^)t 2 ) c = ^S^(°)- ( 16 ) 

Q 

Given the above solution, the experimental properties at any temperature, for example the 
structure factor S q and the dielectric susceptibility as a function of frequency are calculated 
from Xqi 00 ) given by Eq. (llljl . in which A is self-consistently determined by Eq. ( 1T41 and the 
static polarization ((ui)^) self-consistently through Eqs. ATS]) and (ITBjl . 

We now present results of the calculations to illustrate the physical principles as well 
as compare with experiments. Here, the model parameters are obtained from fits to the 
experimental structure factor of PMN (see Fig.H]). Fig. [1] presents the calculated temperature 
dependence of the (mean) transverse optical model frequency at q = suitably normalized to 
show its different regimes. The frequency remains finite at all temperatures, thus excluding 
long-range order in the model which includes dipolar interactions and even weak disorder 



(long-range order would be present for the model with short-range interactions alone js]). 
Associated with this is the dielectric response Xo( w ) ; shown in Fig. El where T cw and T B 



are also identified (see inset in Fig. [2]). Figs. |3] (a)-(b) give the variation of Tew, Tb, and 
T max with the parameters of the model. We see that the characteristic temperatures are 
understood by the ratio of the disorder distribution to the transverse dipole interaction with 
the damping T playing a major role for T max . 

We now compare the static structure factor S q with that measured by neutron scattering 
experiments [sllisj]. S q is derivable from Xgi 00 ) D Y the standard procedure. We obtain the 
following result, 



Without compositional disorder (A = 0), S q= o diverges as the temperature approaches the 
critical point T c where long-range ferroelectric order sets in. For finite compositional disor- 
der (A > 0), S q= Q remains finite at all temperatures and there is no long-range ferroelectric 
order. We find this temperature behavior is in agreement with that observed by neutrons 
in PMN, as shown in Fig. HJ The flat behavior of S q at low temperatures (T — > 0) is due 
to zero-point fluctuations: in the classical limit (H — > 0), S q of Eq. (ITT)) increases with 
decreasing temperature and is finite at T = 0. 

Fig. H] (b)-(c) compare the wavevector distribution of S q with that observed for the 
relaxor PMN at T = 150 K and T = 300 K. The observed line shape cannot be described by a 
simple Lorentzian seen for conventional perovskite ferroelectrics or the "squared Lorentzian" 
expected for random field models. In a supplementary section, we give the details of the 
spatial correlation functions obtained from S q including their anisotropy, which may be 
further tested in experiments. We also discuss the concept of polar nanoregions, which is 
often invoked to qualitatively explain the properties of relaxors. 

Finally, we summarize the conclusions of this work. We have used the insight that the 
frustrating nature of the dipolar interactions (due to their anisotropy) introduced in a per- 
ovskite due to a putative displacive transition together with quenched disorder impedes a 
ferroelectric transition at all temperatures and leads instead to a region of extended ferro- 
electric fluctuations. The nature of the dipolar interactions is such that the physics cannot 
be captured in a mean-field like or two body correlation approximations to the problem. 
Using a minimal approximation scheme which is able to handle the special nature of the 
problem, we are able to derive the observed structure factor of the relaxors and relate it to 
their static and dynamic microscopic properties. 




(17) 
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FIG. 1. Temperature dependence of the zone center optic phonon frequencies Qq for dipole inter- 
actions and compositional disorder. The frequencies remain finite down to T = for finite disorder, 
thus excluding long-range ferroelectric order. 
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FIG. 2. Temperature and frequency dependence of the real part of the dynamic susceptibility xo( w )- 
Inset: Inverse of the real part of the static susceptibility with compositional disorder (A 2 /(t>jj~a) 2 = 
0.03). Deviations from Curie- Weiss law are indicated by the gray dashed line. Blue, red, and green 
lines correspond to ooT/v^ = 2.0 x 10" 2 , 3.0 x 10~ 2 , 4.0 x 10~ 2 . 




9 




FIG. 3. (a) Curie- Weiss (Tew) an d Burns (Tb) temperatures dependence on compositional dis- 
order, (b) Frequency and disorder dependence of the temperature T max at which the dielectric 
constant is maximum. 




FIG. 4. (a) Temperature and (b)-(c) wavevector dependence of the structure factor S a . Red 

nn 

dots and blue squares correspond to neutron scattering data from Refs. [5J, |13[ in the vicinity of 
the (110) and (100) Bragg peaks in PMN. Black solid line corresponds to the calculated S q with 
(vfr-K)/vfr = 0.85, ja 2 /v^ = 0.09, h/(Mv^a A ) l l 2 = 5.0, A 2 /(v^a) 2 = 0.03, and k B T c /(v^a 2 ) = 
1.0 with T c = 465 K. Green dashed line is a Gaussian Bragg peak of width given by the experimental 
resolution used in Ref. [fj (~ 0.01 rlu). 
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In this supplement, we compute the spatial dependence and anisotropy of the correla- 
tion functions of polarization at various characteristic temperatures and various normalized 
disorder strengths. These spatial correlations have been Fourier transformed to give the 
static structure factor S q which has been described in the text and which we have shown 
compares well with the experimental results from neutron scattering. One of the purposes of 
this section is to contrast the correlations of our model to those expected from hypothetical 
polar nanoregions (PNRs) on which we will comment at the end. 

The correlation functions of polarization, {(uiUj) T / , are obtained by the Fourier trans- 
form of the structure factor, 



Sq = ^—coth(^)+ - 2 . (1) 

(MQ q ) 



2MQ q 



We first consider the correlation functions of polarization at large distances (|-Rjj| — > oo). 
For low temperatures (T — y 0) and arbitrarily small but finite disorder (0 < A/ (vqCl) <C 1), 
the correlations are anisotropic and slowly decaying functions, 



, (T ) 2 7^Y x if {ftp \ R ^\ -+ °°> ^ II * (longitudinal), 



) 2 7^ x \ — |.Ri 7 -| -> oo, Rn ± z (transverse) 

v. 3 / (v-i-Y VBZ I-^jI 



where (xo ) 1 = (xo) 1 /vq is a. small but finite dimensionless inverse susceptibility; 

and vbz is the volume of the Brillouin zone. The ratio of the longitudinal to transverse 

Q In 

components is proportional (in absolute value) to (xq) , indicating that the positive lon- 
gitudinal correlations are stronger than the negative transverse components. This behavior 
is similar to that of ferroelectrics without random fields. For an uniaxial ferroelectric with- 
out compositional disorder and above the critical temperature T c , the correlation functions 
exhibit the same anisotropy, power law decay, and longitudinal to transverse ratio The 
correlations of Eq. (j2J) are, however, in sharp contrast with those of short-range interactions 
with quenched random fields in three dimensions where the correlation functions decay ex- 
ponentially j^]. 

We now discuss the near-neighbor correlations with compositional disorder. Here, we 
use the model parameters given in Fig. 4 of the paper. For temperatures slightly above T c , 
the near-neighbor longitudinal correlations with and without disorder show similar decay 
with an overall strength that decreases with increasing compositional disorder (Fig. Ufa)). 



This behavior persists down to T = for finite compositional disorder (Fig. QJb)). For no 
disorder, the correlations exhibit the long-range order expected at T = 0. The correlations 
in the transverse direction follow similar decay as that of the longitudinal components except 
that they are negative and significantly weaker (Figs. [Tfc)-(d)). 

To compare with the idea of PNRs invoked to rationalize the behavior of relaxors, we note 
that we find that significant correlations develop below the Burns temperature Tb- However, 
we do not find a difference in the characteristic behavior at short distances compared to that 
at long distances, which might have been expected for PNRs. We find that the correlations 
have a power law behavior which joins smoothly to a short-range part where they must 
saturate to near the on-site correlations. This is characteristic of the fluctuation regime of 
any cooperative problem. We have emphasized that relaxors may be looked on as materials 
in which the fluctuation regime extends from the Burns temperature all the way to T = 0. 

We cannot therefore support the qualitative picture of PNRs. They are usually invoked in 
a suggestive fashion and used to obtain some but not all the principal properties of relaxors 
in Debye-type model of relaxing entities Since we are able to explain all such properties 
as well as S q from our model and calculations, we feel on secure grounds. 
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FIG. 1. Log-linear plot of near-neighbor correlation functions of polarization \(uiUj) th ) for 
several temperatures and disorder strengths. 
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